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Abstract 

We study the possibility of observing the stau signal at LHC in case that the gravitino 
is the lightest supersymmetric particle and the stau is the next lightest supersymmetric 
particle in high reheating temperature scenario. We show that a number of stau signals 
can be observed at LHC for y/s = 7 TeV and an integrated luminosity of 1 fb~^ in most 
of the parameter region for the reheating temperature T]i>O{10^) GeV. We also show 
that the parameter region with Tfl>2 x 10^ GeV, which is consistent with the thermal 
leptogenesis, is all covered for = 14 TeV with 10 fb~^. 



1 Introduction 



Supersymmetry (SUSY) is one of the most plausible models beyond the standard model. In 
some classes of minimal supersymmetric standard model (MSSM) with the gravitino, the 
gravitino (G) becomes the lightest supersymmetric particle (LSP) and the lightest stau (f ) 
becomes the next-to-lightest supersymmetric particle (NLSP). Such a scenario may predict 
the stau lifetime longer than 0{1) sec because of the weak coupling between the gravitino and 
the stau. The existence of the long-lived charged massive particles like the stau is appealing 
from the viewpoint of the discovery at the collider, since they are observed as charged tracks 
at the detector. The search for such long-lived particles at LHC has already started [1]. 

However, there are some cosmological problems in this scenario. The present energy 
density of the gravitino produced in the thermal scattering at the epoch of reheating may 
exceed the observed dark matter energy density. Moreover, the existence or the late-time 
decay of the stau at the epoch of Big Bang Nucleosynthesis (BBN) may spoil the successful 
BBN. It has been pointed out that, from these two cosmological constraints, upper bounds on 
the gluino mass are obtained for a given reheating temperature and for a given stau mass [2] 
(see also [3]). Such upper bounds on the gluino mass predicts promising collider signatures 
of this long-lived stau scenario. 

In this paper, we study the observability of the stau signals at LHC in this long-lived 
stau scenario under the cosmological constraints, assuming that the R-parity is conserved, 
and that there is no entropy production after the reheating. We show that a higher reheating 
temperature predicts a lower gluino mass, and therefore more stau signals can be observed 
at LHC. It is shown that a number of stau signals can be observed at the early LHC (the 
center-of-mass energy y/s = 7 TeV and the integrated luminosity Lj = 1 fb^^) in most of 
the parameter region for Tr>C'(10^) GeV. We also show that the parameter region with 
Tr>2 X 10^ GeV, which is consistent with the thermal leptogenesis [4, 5], is all covered for 
^/s = 14 TeV and Li = 10 fb~^. This paper is the complete version of our previous work [6]. 
In addition to the analysis in our previous paper, we include the followings in this paper. We 
study the parameter region where the stau annihilates near the pole of CP-even heavy Higgs 
boson. The upper bounds on the gluino mass are shown in various gaugino mass relations. 
The results using the detector simulation are shown in more details. 

This paper is organized as follows. In Section 2, we discuss the cosmological constraints. 
In Section 3, the constraints from Tevatron and the signatures at LHC are studied. We 
conclude this paper in Section 4. 
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2 Cosmological Constraints 



In this section, we discuss the cosmological constraints in the scenario with a long lived stau 
and with high reheating temperature. There are mainly two constraints: the constraints 
from gravitino over-production and the constraints from BBN. The over-production bound 
is reviewed in Sec.2.1, where an upper bound on the gluino mass is obtained for a given 
gravitino mass and a reheating temperature. In Sec.2.2, we discuss the BBN bound, and 
we get an upper bound on the gravitino mass for a given stau mass. Combining these two 
cosmological constraints, an upper bound on the gluino mass is obtained for a given stau 
mass and a reheating temperature. 

We discuss the three different annihilation processes of staus: (i) electrowcak processes, 
(ii) annihilation into light Higgs bosons(/i) [7, 8] and (iii) annihilation near the pole of CP-even 
heavy Higgs boson(if) [8]. 

2.1 Gravitino Over-production Bound 

The gravitinos are produced by the scattering process of particles in thermal bath after the 
epoch of reheating [9, 10, 11, 12]. The gravitino abundance takes the form of [11] 
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where 7713/2, rn~, and m~ are the physical masses of the gravitino, the bino, the wino 
and the gluino respectively. The reheating temperature Tr is defined by 



where is the inflaton decay rate, is effective degrees of freedom and Mp = 2.4 x 10^^ 
GeV is the Planck scale. We used the one-loop renormalization group equations to evolve 
the running masses of gauginos up to the scale jj, = Tr. The numerical coefficients, which 
depend on Tr logarithmically, are evaluated at Tr = lO^GeV in Eq.(l). In the numerical 
analysis, we include those logarithmic dependences. Note that Eq.(l) potentially includes an 
0{1) uncertainty [10, 12]. 

Although gravitinos can be also produced by inflaton decay [13] and the moduli decay [14], 
they are model dependent, and we do not include these production processes, for simplicity. 
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In the present scenario, late-time stau decay also produces gravitinos, 

or,?*^.2.8xlO-(Ji,)(ig^). (3) 

where Y~ = n~/s is the stau abundance, n~=n~ +n~ is stau number density after the freeze- 
out, and s is the entropy density. However, its contribution is negligible in the parameter 
region of our interest, and we neglect it in the following discussion. 

The gravitino energy density should not exceed the observed DM density [15], 

^i/2h^ < 0.122 (95%C.L.). (4) 

This constraint gives the upper bound on the value of gluino mass for a given gravitino 
mass and reheating temperature, for fixed values of bino mass and wino mass. The most 
conservative bound is obtained for = = m3/2, which is shown in Fig.l. In the 
lower gravitino mass region, the term proportional to m^/m^i2 dominantly contributes to 
the right-hand side in Eq.(l). Thus, in that region, the upper bound on m~ is proportional to 
^m3/2, as we can see in Fig. 1. On the other hand, we see that the upper bounds drop down 
in the higher gravitino mass region, because the terms proportional to m3/2 and m~^^/m3/2 
become non-negligible. 

We find that in gravitino LSP scenario, the gluino mass cannot be higher than 2.5 TeV 
for Tji>2x 10^ GeV which is required in the thermal leptogenesis [4, 5]. It is also found 
that the reheating temperature higher than 0(10^°) GeV is not allowed in the scenario 
with gravitino LSP. Note that these bounds are the most conservative ones and should be 
satisfied independently of the BBN constraints, unless there is an entropy production after 
the gravitino production. 

2.2 BBN Bound 

In the present scenario, staus are long-lived so that they may affect BBN. There is an upper 
bound on the stau lifetime depending on the value of stau relic abundance at the BBN epoch. 
The stau lifetime is given by 

f~= ^JZ^ , (5) 

and hence the upper bound on the stau lifetime can be translated into an upper bound on 
the gravitino mass for a given stau mass. 
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Figure 1: The upper bound on the gluino mass in case that M~ = = 771,3/2 • 

Among the constraints on stau hfetime from BBN, the most stringent bound comes from 
the catalyzed effect where stau forms a bound state with ^He and overproduces ®Li [16] . The 
constraint is roughly r~ < 1000 sec for Y~> 10^^^, while there is almost no constraint on stau 
lifetime for 10^^^. It is also important for large Y~ that the energetic hadrons produced 
by stau-decay modifies the abundance of D nuclei [17]. 

Stau abundance Y~ depends on the value of the stau annihilation cross section, and hence 
on the stau annihilation process. In this paper, we consider the following three different 
parameter regions. 

(A) In most of the parameter region of the MSSM, the stau annihilation is dominated by 
the electroweak interaction. 

(B) When the f — f — h coupling is large, the stau annihilation into the light CP-even Higgs 
bosons h is enhanced [7, 8]. 

(C) When the f — f — H coupling is large and the heavy Higgs mass satisfies the condition 
rriH ~ 2m~, the staus can annihilate at the resonance of the heavy Higgs boson H [8] . 

In the cases of (B) and (C), the stau abundance Y~ can be significantly reduced compared 
to the case (A) [7, 8]. In the following subsections 2.2.1-2.2.3, we discuss the gravitino mass 
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upper bound for a given stau mass in these three different parameter regions. The obtained 
BBN bound on the gravitino mass is combined with the over-production bound discussed in 
Sec.2.1, which leads to the gluino mass upper bound for a given stau mass and reheating 
temperature. 

2.2.1 (A) Stau Annihilation via Electroweak process 

When the electroweak process is dominant, the stau abundance is estimated as [18] 



where the constraint in the region lOOGcV < m~ < 450GcV comes from the bound on Li 
overproduction by the catalyzed effect, whereas the bound on the hadronic decay of staus 
gives the constraint in the region 450GeV < m~ < lOOOGeV. 

This upper bound on the gravitino mass can be easily translated into the upper bound 
on the gluino mass for a given stau mass and a reheating temperature. As we can see from 
Eq.(l) and Eq.(4), the upper bound on the gluino mass also depends on the other gaugino 
masses. In Fig. 2, we show the upper bound on the gluino mass for a given gaugino masses 
and reheating temperature. We set the stau mass to be 300 GeV in the figure. We find the 
heavier gaugino masses lead to the more stringent upper bound on the gluino mass. 

In the following discussion, we show the upper bounds on the gluino mass for the following 
two different gaugino mass relations: 

(i) = = l.lm~. This gives almost the most conservative upper bound on the gluino 
mass (cf. Eq. (1)). Note that if is too degenerate with m~, the stau abundance 
is enhanced by the B{W) decay, resulting in more stringent bounds on Tr (or m~) [18]. 

(ii) Gaugino masses satisfy the GUT relation, rn-/ gl{m~) = m^/g2{m^) = Tn~/g'^{m~), 
where (/x) , (a*) and 53(/u) are the running gauge coupling constants of U{1)y, SU{2) l 
and SU (3)c gauge symmetries at the mass scale n respectively. This is the case for the 
minimal supergravity and the minimal gauge mediation models. 

Fig. 3 and 4 show the results for the cases (i) and (ii), respectively. Prom the figures, we find 
the foUowings; 




(6) 



Then, the bound from BBN on gravitino mass is [17] 



^ f 0.4 GeV - 10 GeV (lOOGeV < m~ < 450GeV) 
'3/2 ~ I 10 GeV - 20 GeV (450GeV < m~ < lOOOGeV) 



(7) 
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Figure 2: The upper bound on the gluino mass for the case (A), where the stau annihilation 
is dominated by the electroweak processes. The bino and wino masses are varied assuming 
m~ = m^, while the stau mass is fixed as m~ = SOOGeV. The solid (blue) lines are upper 
bounds on the gluino mass. We also plot the number of stau signatures at LHC with each ^/s 
and integrated luminosity. Here, the events only include the productions of gluinos, charginos 
and/or neutralinos (and staus), and we impose cuts and triggers. The horizontal dashed line 
around m~ ~ 470 GeV comes from the CDF bound. We discuss the collider signature in 
Section. 3. 

(i) m~ = = l.lm- case: the gluino mass cannot exceed about 2.4 TeV for m~ < 
1 TeV and rR>10'^ GeV. For Tr > 3 (5) x 10*^ GeV, the stau mass is bounded as 
m~<700 (500) GeV, and the gluino mass is bounded as m~<1100 (700) GeV. There 
is no region for Tr > 7 x 10^ GeV. 

(ii) The case of the GUT relation: the constraint on gluino mass and the reheating temper- 
ature are much severer than the case (i). There is no region for Tr > 7 x 10^ GeV. 



2.2.2 (B) Stau Annihilation with large stau-stau-light Higgs coupling 

When the t — t — h coupling is large, the stau annihilation into the light CP-even Higgs 
bosons h is enhanced [7, 8]. The trilinear coupling of the lighter stau and h is given by 

C = -A-^oTlnh, (8) 
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Figure 3: Case (A)-(i), where the stau annihilation is dominated by the electroweak inter- 
actions, and = mg, = 1.1 m~. The sohd (blue) lines are upper bounds on the gluino 
mass for various reheating temperatures Tr. Contour plots of the number of expected SUSY 
events at LHC and the line of the CDF bound are shown in the same way as Fig. 2 . The 
stau is not the LSP under the dotted (purple) line. 

where the coefficient takes the form of 



-gzMz 



-- + sin 6w cos 6r — sin Owsin Or 



(9) 



where g, gz, M\y, Mz, rrir and 9w are the Standard Model parameters, ^ is the Higgsino 
mass parameter, tan/3 is the ratio of VEVs of the two Higgs doublets, and 9r is the mixing 
angle of the staus defined by 




cos^T- sin^T- 
— sin^T- cos^T- 




(10) 



This trilinear coupling A~f^o is enhanced when fi, tan/3 and sin20T- are large. The stau anni- 
hilation process via the large trilinear coupling may significantly reduce the stau abundance 
[7, 8]. 
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Figure 4: Case (A)-(ii), where the stau annihilation is dominated by the electroweak inter- 
actions, and the gaugino masses satisfy the GUT relation. The contour lines are the same as 
Fig. 3. 

However, the large trilinear coupling may create disastrous charge/color breaking (CCB) 
minima in the t — h potential [7]. In such a case, the vacuum in our universe is a local 
minimum, and therefore, our vacuum will eventually decay into a global minimum. The 
condition that the lifetime of our vacuum must be longer than the age of the universe, gives 
the upper bound on the value of A—f^o, as discussed in Appendix A. The resultant upper 
bound on A~^o is shown in Fig. 5 for nih = 120 GeV. 

The upper bound on A—f^o gives the upper bound on the stau annihilation cross section. 
Since the relic abundance of the stau is given by 



the upper bound on the stau annihilation cross section gives the lower bound on Y~, where a 
is the annihilation cross section of stau, v is the relative velocity of staus and () denotes the 
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Figure 5: CCB bound on the stau-Higgs coupling .4~^o for = 120 GeV. The constraint 
from the zero temperature decay rate is severer for m~ > 220 GeV, while the finite tempera- 
ture transition dominates the bound for m~ < 220 GeV. 



thermal average. In the large trilinear coupling case, {av) takes the form of 

(^-) - ^fe A + ^^itzf ft^ (12) 



where 



yi-rfe (l-n/2Vl-rf 

Jf'^n 7r^^^^~^h), ft = 7z 7TT2 0{l-rt), (13) 

(1 - (1 - Th 4:y 



with rfi = m|/m~ and rt = mf/m~. If m~ is larger than m^, the annihilation into Higgs 
bosons via the large trilinear coupling is possible, and the first term in Eq.(12) becomes 
nonzero. In the same way, the second term in Eq.(12) becomes nonzero when m~ is larger 
than mt- In the following, we take m/j = 120 GeV and rrit = 173 GeV. The lower bound on 
in this large A~^o scenario is shown in Fig. 6. For comparison, we also show the Y~ when the 
electroweak process is the dominant stau annihilation process in the same figure. We find that 
Y~ in the large -4~^o scenario is significantly smaller than that of the electroweak annihilation, 
but it cannot be smaller than 10~^^. Therefore, the bound from the ^Li overproduction by 
the catalyzed effect severely constrains the gravitino mass [17], 

^3/2 < 0.4GeV - 115GeV (lOOGeV < m~ < lOOOGeV) (14) 

as shown in Fig. 7. 
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Figure 6: Green solid line: Y~ vs m~ when A~f^ is the maximal value. We take = 120 GeV 
and mt = 173 GeV. Blue dashed line: Y~ when the electroweak process is the dominant stau 
annihilation process. 



This gravitino mass upper bound leads to the upper bound on the gluino mass for a given 
stau mass and reheating temperature. As in the case in Section 2.2.1, we show the upper 
bounds on the gluino mass with the two gaugino mass relations:(i) m~ = niy^ = l.lm~, and 
(ii) the GUT relation. The upper bound with the relation (i) is shown in Fig. 8, while that 
with the relation (ii) is shown in Fig. 9. From Fig. 8 and Fig. 9, we find the followings; 

(i) m~ = = l.lm~ case: the behavior of the gluino mass upper bound in Fig. 8 are 
similar to that in Fig. 3 for m~ < 450 GeV, since the gravitino mass upper bound is 
determined by the constraints from ^Li in the region. On the other hand, the behavior 
is quite different for m~ > 450 GeV, because the gravitino mass upper bound in the 
case (A) is determined by the constraints from the deuterium, while the upper bound 
in the case (B) is still determined by the constraints from ^Li. It is also remarkable that 
Tr > 10^ GeV is possible. 

(ii) The case of the GUT relation: the behavior of the gluino mass upper bound is almost the 
same as Fig. 4. We also see that the behavior of the gluino mass bound slightly changes 
at m~ = nit, because the annihilation into top quark becomes possible for m~ > mt- 
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Figure 7: Green solid line: gravitino mass bound for the case with an enhanced /i-stau 
coupling . Blue dashed line: gravitino mass bound for the case of the electroweak annihilation. 

2.2.3 (C) Stau Annihilation near the pole of the Heavy Higgs Boson 

The stau annihilation cross section can be considerably enhanced in case that niH — 2m~ 
and At, tan/3 and sm29r are large [8]. The trilinear coupling of the lightest stau with H 
takes the form of 

C = -A~HnnH, (15) 

where the coefficient is given by 

in the large tan/3 region. When Arianji is large, the first term in the right hand side becomes 
enhanced. In such large Artanfi case, the annihilation through H exchange can be the 
dominant process, and if mn satisfies the resonance condition niH — 2m~, it is enhanced 
greatly. 

The large -4—^ may generate the disastrous CCB minimum, and therefore, the value of 
.A—ji^ is constrained as in the case of the enhanced /i-stau trilinear coupling. We obtain the 
upper bound on A~^ as we discuss in Appendix A. The bound on A—^ is shown in Fig. 10 
for rriH = 2m~. 

In this case, the stau relic abundance can be sufficiently reduced to avoid the BEN con- 
straint, i.e., Y~ < 10^^^. For illustration, we take niH = 2m~ and ^—^^^=1.3 x m~, which 
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Figure 8: Case (B)-(i), where the stau annihilation is enhanced by the /i-stau-couphng, and 
= = 1.1 m~. The contour hues are the same as Fig. 3. 

is below the upper bound on A~jj in lOOGeV < m~ < lOOOGeV, as can be seen from 
Fig. 10. The resultant abundance Y~ is shown in Fig. 11, where we use the program mi- 
crOMEGA2.4 [22] to calculate the relic abundance. We see that Y~ can be less than 10""^^ 
for lOOGeV < m~ < lOOOGeV, and therefore BBN does not provide any constraints on the 
gravitino mass. 

There is another constraint from the CMB spectrum distortion [23]. However, the dif- 
ference between the stau mass and the gravitino mass upper bound m^2^ is small, (m~ — 
m^2^(m~))/m~ < 0.01. Thus, we take my'^{m~) = m~, for simplicity. 

The upper bound on the gravitino mass leads to the upper bound on the gluino mass for 
a given stau mass and reheating temperature. As in the case in Section 2.2.1 and Section 
2.2.2, we show the gluino mass upper bound (i) for m~ = = 1.1 x m~ in Fig. 12, and 
(ii) when there is the GUT relation between gaugino masses in Fig. 13. Note that, in the 
larger stau mass region, the gluino mass is maximized when the gravitino mass is smaller 
than m^2^'(m~). That is because the upper bound on the gluino mass drops down in the 
higher gravitino mass region as shown in Fig.l. From Fig. 12 and Fig. 13, we find followings; 
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Figure 9: Case (B)-(ii), where the stau annihilation is enhanced by the /i-stau-couphng, and 
the gaugino masses satisfy the GUT relation. The contour lines are the same as Fig. 3. 

(i) m~ = = 1.1 X m~ case: Interestingly, Tr > 2 x 10^ GeV which is required by 
thermal leptogenesis is possible in the region m~<1500 GeV. We can even achieve 
Tr > 5 X 10^ GeV. The gluino mass bound becomes lower in the larger m~ region 
because the Bino mass and the Wino mass becomes heavier for large stau mass, and the 
constraints of the gravitino overclosure becomes severer. 

(ii) When there is the GUT relation: Tr > 2 x 10^ GeV is possible in m~<350 GeV. The 
gluino mass is bounded as m~< 2 TeV when Tr > 2 x 10^ GeV. 

CD Q -I L 

3 Collider Signatures 

Let us study the collider detectability of the scenario of the high reheating temperature. In the 
previous section, it was shown that the gluino mass is constrained to be less than a few TeV 
to realize a high reheating temperature. Light SUSY particles have been already excluded by 
direct searches in Tevatron, while TeV colored SUSY particles are in the reach of sensitivity 
of LHC. If heavy SUSY particles are produced at a collision, they subsequently decays into 
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Figure 10: CCB bound on the stau-Higgs coupling The finite temperature transition 

dominates the bound for m~ > 100 GeV. 

hghter SUSY particles with radiating SM particles, and promptly generate the NLSP stau 
in the end of the decay chain. Since the gravitino mass is as large as 0(1 — 100) GeV in the 
high reheating temperature scenario, the NLSP stau is long-lived enough to be observed as 
a stable particle in the detectors. Noting that the stau has an electromagnetic charge, we 
expect the events with charged tracks when SUSY particles are produced at collisions. In 
this section, we discuss the Tevatron bound and the LHC sensitivity of the high reheating 
temperature scenario. 

Before proceeding to the collider study, we summarize the tools for the numerical analysis. 
We use PYTHIA 6.4.22 [24] to study the kinematics and to estimate the cross sections except 
for those of the colored SUSY particle productions. The gluino and squark production cross 
sections are estimated by the program Prospino2 [25] at the NLO level. The Tevatron and 
LHC detectors are simulated by the package PGS4 [26]. 

3.1 Tevatron Bound 

At the Tevatron experiments, the stau is expected to behave as a heavy muon, namely a 
charged massive and long-lived particle. As long as the stau has a large velocity, it is not 
distinguishable from the muon. When the transverse momentum, px, of the muon is large, 
it has a large velocity of /3 ~ 1, while the stau is likely to have a lower speed. Such a high-pT^ 
and low speed "muon" has been searched for in Tevatron by measuring the time of flight 




15 



1e-15r 

9e-16 

8e-16 

7e-16 

6e-16 




100 200 300 400 500 600 700 800 900 1000 

rur (GeV) 
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[27, 28]. According to [27], the events are selected by the following trigger: 

• the highest px "muon" candidate has px larger than 20 GeV which satisfies an isolation 
condition ET{0.4:)/pT{fJ.) < 0.1, 

where £'t(0.4) is the sum of the transverse energy within a cone R = 0.4 around the candi- 
date, excluding the energy deposited by the muon candidate itself, and prifJ-) is the transverse 
momentum of the highest pT muon candidate. Note that the long-lived charged massive par- 
ticle is identified as a "muon" candidate. The events that satisfy the above trigger condition 
are read out. The Tevatron constraint is that the production cross section of the long-lived 
charged massive particle which runs toward the direction \ri\ < 0.7 with pT > 40GeV, and 
with the velocity 0.4 < (3 < 0.9, is smaller than 10 fb at the 95% C.L.. The constraint gives 
the lower bound on masses of SUSY particles. 

We show the constraints from Tevatron in Figs. 2-4, Fig. 8, Fig. 9, Fig. 12 and Fig. 13. The 
constraints in Fig. 3, Fig. 8, and Fig. 12 are the same, which corresponding to the cases (A)- 
(i), (B)-(i), and (C)-(i), respectively. Similarly, the constraints in Fig.4, Fig. 9, and Fig.l3 
(corresponding to the cases (A)-(ii), (B)-(ii), and (C)-(ii), respectively) are the same. We only 
consider productions of the gauginos and the lighter staus, and do not consider the production 
of the other scalar particles. This is realized when their masses are relatively heavy. In the 
case that the gluino mass is relatively small, gluino pair production is the main production 
channel. That is the case in the regions around the horizontal lines in Fig. 2, Fig. 3, Fig. 8 and 
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Figure 12: Case (C)-(i), where mn — 2m~ and the i/-stau-couphng is enhanced, and = 
= 1.1 m~. The contour hues are the same as Fig. 3. 

Fig. 12. When the gluino mass is large, the stau direct production, chargino-neutralino and 
chargino-chargino pair productions are the main production channel. The vertical lines of the 
CDF constraints in Fig. 3 are determined by the chargino-neutralino and chargino-chargino 
pair production cross sections. The same situation holds in the regions around the vertical 
lines in Fig. 8 and Fig. 12. On the other hand, the vertical line in Fig. 4 is determined by the 
stau direct production. That is also the case in the regions around the vertical lines in Fig. 9 
and Fig. 13. 

3.2 LHC signatures 

The scenarios of the high reheating temperature predict the gluino mass being less than a 
few TeV. Although the Tevatron energy is not large enough to cover the mass range, the 
LHC is suited for detecting the particle. The LHC is running at the center-of-mass energy 
^/s = 7 TeV and the integrated luminosity is planned to become up to afewfb^^ in 2011. The 
schedule of 2012 is still in discussion: an optimistic scenario is ^/s = 8 TeV with Li ~ 10fb~^. 
After the upgrade, the collider is aimed to run at ^/s = 14 TeV. In the following, we 
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Figure 13: Case (C)-(ii), where mn — 2m~ and the if-stau-couphng is enhanced, and the 
gaugino masses satisfy the GUT relation. The contour hnes are the same as Fig. 3. 

discuss the LHC sensitivity of the high reheating temperature scenario in the three setups: 
^ = 7 TeV with Li = 1 fb"\ ^/s = 8 TeV with U = 5 fb"^ and ^/s = U TeV with 
Li = 10 fb~^ 

Analogous to the Tevatron, the stau signature is the charged track of the low speed muon- 
like particle with high px within the detectors. First of all, the SUSY events are selected 
by imposing a trigger menu, and then analyzed off-line with cut conditions in order to be 
distinguished from the background. In the following analysis, we impose the following trigger 
conditions, 

• at least one isolated electron has px > 20 GeV, 

• at least one isolated muon has pi > 40 GeV, 

• at least one isolated tau has > 100 GeV, 

• at least one jet has px > 200 GeV, 

• at least three jet has px > 100 GeV, 

• at least one isolated stau has px > 40 GeV within the bunch. 
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• at least two staus have pt > 40 GeV within the bunch. 

If any one of these conditions is satisfied, the event is read out. In our simulation, the isolation 
conditions on the electron and the tau relies on PGS4, while those of the muon and the stau 
are 

1. the summed in a R = 0.4 cone around the particle (excluding the particle itself) is 
less than 5 GeV, 

2. the ratio of Et in a 3 x 3 calorimeter array around the particle (including the particle's 
cell) to pt of the particle is less than 0.1125. 

The bunch condition, i.e., the condition that the stau has a velocity large enough to reach 
the muon trigger detector before the next bunch collides, is necessary for the stau triggers, 
since otherwise they do not work correctly [30]. We require /3 > 0.7 for the stau propagating 
in the barrel region (|?7| < 1.0) and /3 > 0.8 in the endcap region (1.0 < |?7| < 2.8) [29]. 

Among the events that are read out by the triggers, the SUSY signals are required to 
satisfy the following cut conditions, 

• PT > 20 GeV 

• 0.5 < /3 < 0.9 

• \r)\ < 2.5. 

Although the muon productions arc the relevant standard model backgrounds of the stau 
signals, they are significantly reduced to be almost zero by especially the first two cuts, since 
the high pr muons have /? ~ 1. The last condition is added because the stau is detected in 
the muon detector [29, 30]. 

The number of events is reduced to roughly 50% ~ 90% by the above trigger conditions. 
When the stau direct production or chargino/neutralino pair production is the dominant 
production channel, roughly a half of the triggered events are selected by the stau and the 
others are selected by jets and taus. On the other hand, when the gluino production is 
dominant, most of the triggered events are selected by jets and taus. Note that we put the 
trigger efficiency as unity for simplicity. In the real detector system, however, the efficiency 
for single jet and for single tau is not good. Thus the number of the events selected only by 
single jet or single tau (roughly 10% ~ 40% of the triggered events) may be further reduced 
if we consider the trigger efficiency. 

We also have to consider the reconstruction efficiency of the stau. The reconstruction 
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efficiency varies from 0.1 to 0.9 for /3 ~ 0.5 ~ 0.9 according to ATLAS CSC studies [29]. The 
analysis can be improved by the method studied in [30], which provides the efficiency more 
than 90% for /3 > 0.5. In the following, we assume the efficiency for the stau with ^ > 0.5 to 
be 100% and for the stau with ^ < 0.5 to be zero for simplicity. 

In Fig. 2-4, Fig. 8, Fig. 9, Fig. 12 and Fig. 13, we show the number of stau signals that satisfy 
the trigger conditions and the cut conditions at LHC for ^/s = 7 TeV with Lj = 1 fb~^, 

= 8 TeV with Lj = 5 fb^"*^, and ^/s = 14 TeV with Lj = 10 fb~^. The contour lines in 
Fig. 3, Fig. 8, and Fig. 12 are the same, which corresponding to the cases (A)-(i), (B)-(i), and 
(C)-(i), respectively. Likewise, the contour lines in Fig. 4, Fig. 9, and Fig. 13 (corresponding to 
the cases (A)-(ii), (B)-(ii), and (C)-(ii), respectively) arc identical. We only consider gaugino 
production and stau direct production at the collision as the production channels, and do 
not consider the production of the other scalar particles. As in the case of the Tevatron 
constraints, the gluino pair production is the main channel if the gluino mass is relatively 
small. That is the case in the regions around the horizontal lines in Fig. 2, Fig. 3, Fig. 8 and 
Fig. 12. When the gluino mass is large, the stau direct production, chargino-neutralino and 
chargino-chargino pair productions are the main production channel. The vertical lines of 
the LHC signatures in Fig. 2 arc determined by the chargino-neutralino and chargino-chargino 
pair production cross sections. The same situation holds in the regions around the vertical 
lines in Fig. 3, Fig. 8, Fig. 12, and the horizontal lines in Fig. 4, Fig. 9, Fig. 13. On the other 
hand, the vertical line in Fig. 4 is determined by the stau direct production. That is also the 
case in the regions around the vertical lines in Fig. 9 and Fig. 13. 

Now let us discuss the implications of these LHC signatures on the reheating temperature. 
Combining the cosmological upper bound on the gluino mass discussed in Sec. 2 with the 
results obtained in this section, we can see the minimal number of stau signals in each 
reheating temperature. For example, if Tr is larger than 3 x 10^ GeV and the main stau 
annihilation process is the electroweak process, more than 10 staus are observed for -s/s = 
7 TeV with Lj = 1 fb~^. Also, more than 10 staus are observed for ^/s = 14 TeV and 
Li = 10 fb~^, if Tr is larger than 10^ GeV, m~ < 1000 GeV, and the electroweak process is 
dominant. In TABLE 1, we summarize the range of the reheating temperature with which 
more than 10 stau signals are expected in each scenario, y's, and Lj. It is remarkable that 
Tr > 2 X 10^ GeV which is required by the thermal leptogenesis is all covered at s/s = 14TeV 
and Li = 10 fb-^ 

Although we have considered only gaugino production and stau direct production at 
the collisions when discussing the Tevatron constraint and LHC signatures, gluino-squark 
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Table 1: The range of the reheating temperature with which more than 10 stau signals are 
expected in each scenario. 

production becomes dominant when the squark mass is relatively light. We show the Tevatron 
constraint and LHC signatures in Fig. 14 when the squark mass m~ equals to m~ and the 
gaugino masses satisfy the condition (i), m~ = m^ = l.lm~. We also show the upper bound 
on the gluino mass in the case of (A)-(i). The cross section of gluino-squark pair production 
is roughly ten times larger than the gluino pair production. Thus, in the region where the 
gluino pair production is dominant in Fig.3, the expected number of stau signals increase in 
this light squark case. As can be seen from Fig. 14, the contour lines of the constant stau 
number shift upwards by Am- 2± 200 ~ 400 GeV, compared with Fig. 3. On the other hand, 
the result does not change very much when there is the GUT relation even for m~ 2± m~. 
That is because the gluino production is not the dominant channel of the SUSY events in 
the regions around the vertical lines and the horizontal lines of, e.g., Fig.4. 

4 Conclusion and Discussion 

In this paper, the observability of the stau signals at LHC has been studied in the scenario with 
the gravitino LSP and the stau NLSP under the cosmological constraints. It was seen that a 
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Figure 14: Same as Fig. 3 but with mgquark — ^3- 



higher reheating temperature predicts a lower gluino mass, and therefore more stau signals 
can be observed at LHC. The upper bound on the gluino mass depends on the annihilation 
process of the stau, and we have considered three cases that (A) the annihilation is dominated 
by the electroweak processes, (B) the stau annihilation into the light Higgs boson is enhanced, 
and (C) the stau annihilation takes place near the pole of the heavy Higgs. In the case (A), 
which is true in most of the MSSM parameter region, the reheating temperature cannot be 
larger than 7 x 10® GeV. If we further assume that the gaugino masses satisfy the GUT 
relation, the upper bound becomes severer and it cannot exceed 7 x 10'' GeV. On the other 
hand, the reheating temperature can be much higher for the cases (B) and (C). Especially, we 
saw that the reheating temperature can be as large as Tr > 2 x 10^ GeV, which is required by 
the thermal leptogenesis, in the case (C) even if the gaugino masses satisfy the GUT relation. 

We have then investigated the Tevatron constraints and the LHC signatures of the long- 
lived staus by imposing the trigger and the cut conditions. We found that, when the stau 
annihilates mainly via the electroweak process, it is expected that more than 10 stau signals 
are observed at the first stage of LHC with = 7TeV and Li = Ifb^^ in the region 
where Tr is larger than 3 x 10® GeV. When the heavy Higgs boson effectively contributes 
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to the stau annihilation, the first stage of LHC has a sensitivity to the parameter region of 
Tr > 5 X 10^ GeV. The LHC sensitivity improves greatly as the collider energy and the 
luminosity increase. The stau signals will be detected if Tr exceeds 10^ GeV as long as 
m~< lOOOGeV. It is also emphasized that all the regions which satisfy Tr > 2 x 10^ GeV 
can be checked by LHC with ^ = 14TeV and Lj = 10fb~^ 

We have also studied a specific case of the gaugino mass spectrum, i.e. the GUT relation. 
The LHC sensitivity as well as the cosmological constraint on the gluino mass is sensitive 
to the gaugino mass spectrum. The cosmological upper bound becomes severer than the 
case with generic gaugino masses, while at LHC the chargino and/or ncutralino channels 
dominate the productions of the SUSY events instead of the channels of the colored (gluino) 
SUSY particles, especially for a low gluino mass region. It was found that the LHC has a 
detection sensitivity for lower reheating temperature models. 

In this paper, we have assumed that there is no entropy production after the reheating 
epoch. However, our results can also be applied to the case with an entropy production, by 
replacing the reheating temperature Tr with r|? = Tr/A, where A is the dilution factor of 
the gravitino abundance. Even if the entropy production occurs after the freeze-out of the 
stau and before the BBN, the results of the case (C) (Figs. 12 and 13), where there is no BBN 
constraint, hold for . 

It is usually difficult to probe the reheating epoch of the universe directly. Nonetheless, 
the LHC has a sensitivity to the high reheating temperature models, as we discussed in 
this paper. Once heavy charged tracks will be observed at LHC, the reheating temperature 
receives an upper bound, depending on the SUSY mass spectrum. This enables us to reveal 
features of the early universe. 
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A Vacuum stability constraint on the stau trilinear coupling 



In this appendix, the vacuum constraint on the stau trihnear couphng is discussed, which is 
used in Sec. 2.2.2 and Sec. 2.2.3. The hfetime of our vacuum can be estimated by using the 
"bounce method" [20]. By using the euchdian action S'e[^] with the bounce solution (p, the 
decay rate of the false vacuum per unit volume is estimated by 

r/y~£;^exp[-5E[^]] , (i7) 

where E is the typical energy scale of the potential. Wc put the value of the potential energy 
at the global minimum of the potential as 0. Since the lifetime of our vacuum must be longer 
than the age of the universe, 

(18) 

must be satisfied, where is present Hubble constant. This constraint gives the lower bound 
on the bounce, 

SEW\>m. (19) 

Let us start from the case that the f-f-h trilinear coupling .4™^o is large. In our analysis, 
the euclidian action takes the form of 



/oo 
-oo 



\{diT){dif) + \{dih){dih) + c/(f , h) 



(20) 



where f/(r, h) is the potential of r and h. We calculate the potential at one loop level for 
Higgs potential and at tree level for r potential. Now the constraint (19) gives the upper 
bound on ^~;jO for a given stau mass which is shown in Fig. 5 of Sec. 2.2.2, for nih = 120 GeV. 
When we calculate the bounce, we take straight line path from the false vacuum to the true 
vacuum in (r, h) plane. We checked that the full analysis using two dimensional path changes 
the result from the straight line approximation by only 0(1)%. 

We also have to care about the thermal transition of our vacuum to CCB vacuum in the 
early stage of the universe [19]. The procedure of calculating the thermal transition rate is 
similar to that of zero temperature case. We use free energy F instead of euclidian action Se 
to calculate the bounce solution ^. The thermal transition rate at temperature T per unit 
volume is given by 

T{T)/V ~ T^e-^l^'^l/^. (21) 
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Since the thermal transition must not occur, 

, 1 



^*;^(r(r)/F) « 1 (22) 



must be satisfied, where tj and tf are initial and final time respectively, and H(t) is the 
Hubble constant at each time. In our analysis, the free energy takes the form of 



^[T,h,T] = I 



oo 

3 



^{dif){dif) + ^{dih){dih) + u{f,h) + svth{f,h,T) 



(23) 



where SVth is the thermal potential calculated by thermal field theory. The thermal potential 
include the contribution from top quark and gauge bosons of U(l) and SU(2) gauge symme- 
tries in the Standard Model. From the constraint (22), We get the upper bound on .4™/j0- 
The resultant upper bound on .A:^^o for a given stau mass is also shown in Fig.5. 

The upper bound on the f-f-H trilinear coupling which is discussed in Sec. 2.2.3, 

can be obtained in the similar way. The result is shown in Fig. 10 for niH = "^rn-, where we 
take account of the Higgs potential at tree level, and the contributions from bottom quark, 
tau lepton and the gauge bosons of the U(l) and SU(2) gauge symmetries for the thermal 
potential. 
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